Theory of Layered Iron Oxide on Frustrated Geometry: 
Electric Polarization, Magnetoelectric Effect and Orbital State 
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A layered iron oxide i?Fe204 (R: rare-eartii elements) is an exotic dielectric material with charge-order (CO) 
driven electric polarization and magnetoelectric effect caused by spin-charge coupling. In this paper, a theory 
of electronic structure and dielectric property in i?Fe204 is presented. Charge frustration in paired-triangular 
lattices allows a charge imbalance without inversion symmetry. Spin frustration induces reinforcement of this 
polar CO by a magnetic ordering. We also analyze an orbital model for the Fe ion which does not show a 
conventional long-range order. 

PACS numbers: 75.80.+q, 77.80.-e, 75.10.-b, 72.80.Ga 
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Ferroelectric (FE) material has a spontaneous electric polar- 
ization induced by displacement of ion and electronic cloud. 
A subtle balance of long-range Coulomb interaction, lattice 
potential and electron covalency derives the FE transition and 
controls the transition temperature. Recently revived mag- 
netoelectric (ME) effects and related multiferroic phenomena 
are recognized as spin-driven FE transitions. A key 
issue is a non-collinear spin structure in frustrated magnets 
where lattice displacement without inversion symmetry is in- 
duced by the symmetric and/or antisymmetric exchange stric- 
tions. Because of the observed gigantic ME effects, this type 
of materials has been examined from viewpoint of potential 
application as well as fundamental physics. 

When electronic charges are ordered without inversion 
symmetry, a macroscopic electric polarization appears in a 
crystal. This can be another mechanism for ferroelectric- 
ity. This new class of ferroelectricity, which we call "elec- 
tronic" ferroelectricity, is examined recently in the charge 
ordered manganites fl, Q] and several organic salts |0, Ml. 
Rare-earth iron oxides 7?Fe204 (R: rare-earth elements) fs\ 
of the present interest belong to this class of materials. The 
crystal structure consists of paired Fe-O triangular lattices 
[Fig. 1(a)] and R-O block ones alternately stacked along the 
c axis. Since a nominal valence of Fe ion is 2.5+, an equal 
amount of Fe^+ and Fe^+ coexists in the paired triangular 
lattices. In the electron diffraction experiments, Bragg re- 
flections at (/i/3 h/3 3m + 1/2) appear below 320K(= Tqo) 
in LuFe204. This observation indicates a valence order of 
Fe ion, i.e. a charge order (CO) of the Fe 3d electrons f^. 
Around 7co, a spontaneous electric polarization and dielec- 
tric anomalies turn up |10]. Moreover, around the ferrirnag- 
netic spin ordering (SO) temperature (7so =250K) | Til 12], 
the gigantic ME effects are recently discovered iTot 13 Tl4ll . 

These experiments imply that /?Fe204 is not only a CO 
driven FE, but also a multiferroic material due to strong 
charge-spin coupling. A possibility of the ME effects in 
the "electronic" ferroelectrics ubiquitously exist in transition- 
metal oxides and organic salts. However, a guiding principle 
for searching the new class of multiferroic materials based on 
the electronic mechanism has not been fully examined until 
now. In this Letter, we present a theory of a dielectric mag- 



net 7?Fe204 as a electronic ferroelectric and multiferroic ma- 
terial. We address the following issues: (i) origin of the elec- 
tric polarization and the FE transition, (ii) mechanism of the 
coupling between electric polarization and magnetization, and 
(iii) orbital structure of the Fe ion. Present study shows that 
the novel dielectric properties in RFsxOa arise from interplay 
among the geometrical frustration and the multi-degrees of 
freedom of electron. 

Start with the electronic structure in a single Fe ion in 
/?Fe204. An Fe ion is five-fold coordinate with three O anions 
in the plane and two at apices |8]. We calculate the crystalline- 
field splitting for the 3d orbitals in a FeOj cluster based on 
the point-charge model. Split orbitals are assigned to the irre- 
ducible representation for the D^h group: the d^^2_^i (A') or- 
bital, and two sets of the doubly degenerate orbitals, {ady^ — 
bdxy., ad^ + bd^2_^2} {E") and {adxy + bdy^, ad^2_^2 —bd^j] 
(£■') where +b^ = 1 . The degenerate E' orbitals take the 
lowest energy with a — 0.89, although the energy levels for 
E' and £" are close. Thus, each d orbital is singly occupied 
in Fe^+ (5 = 5/2), and one of the lowest degenerate orbitals 
is doubly occupied in Fe^+ (S = 2), which has the orbital de- 
gree of freedom, as well as spin and charge 1 15]. This degree 
of freedom is treated by the pseudo-spin (PS) operator with 
amplitude of 1/2: T; = j ,sdi^^^£,£,'di^' s where t/,-^, is the 
electron annihilation operator with spin s and orbital ^ at site 
i, and (7 are the Pauli matrices. 

Now we set up the model Hamiltonian to describe the elec- 
tronic structures and dielectric properties. The electrical re- 
sistivity in R¥t20A shows an insulating behavior even above 
Tqo llS]. This result implies that the 3d electrons are almost 
localized at Fe ions, and incoherent charge motion occurs by 
thermal hopping between Fe^+ and Fe^+. In such insulating 
systems, the long-range Coulomb repulsion between charges, 
and the superexchange interaction between spins and orbitals 
play major roles on dielectric and magnetic properties. The 
Coulomb interaction Hamiltonian ,3^ — ^ij^^i'^j with the 
electron number n, is mapped onto the Ising Hamiltonian 
= L(i7) ^i;6;'2;- Here, we introduce the PS operator 2/ 
for charge which takes 1/2 (—1/2) for Fe^+ (Fe^+). The su- 
perexchange interaction arises from the virtual electron hop- 
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FIG. 1: (a): A pair of the Fe-O triangular lattices. Thick arro 
represent the Coulomb interactions Va^NN, Vcnn and VcNNN- (b): 
schematic view of CO 1/3 accompanied by the electric polarizatioi 



ping between Fe ions. The Hamiltonian is derived from t 
generalized pd model where the on-site Coulomb and e 
change interactions for the Fe 3d and O 2p electrons and t 
electron hopping fpj between the in-plane NN sites are i 
troduced. Since the Coulomb interactions are larger than t, 
the projection-perturbation theory up to 0{t'^^) is applied {hj 
The high-spin states at Fe ions are chosen as initial states 
the exchange processes, and all possible intermediate stal 
are taken into account. The obtained Hamiltonian is classifi 
by the valences of the NN Fe ions as = J^2 + + 
where J^,m is for the interaction between Fe"+ and Fe" 

Each term includes the several exchange processes denot 

(r) 

by F as J^nm — Y.r '^nm ■ We explicitly show one domin; 
term in ,y^22. 
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Here, I,- is the spin operator with an amplitude of 2, and J^^ 
is the exchange constant. We introduce the orbital PS oper- 
ators Tt[ for the three-kinds of the in-plane nearest neighbor 
(NN) Fe-O bond directions, labeled by / = {a,b,c). These 
are defined by n- — cos(2;rn//3)7j' + sin(2;rn//3)7j-'^ with 
{na,ni,,nc) — (1,2,3). Details in Jtj will be presented else- 
where. It is clearly shown that spin I,-, charge Q'j and orbital 
tt/ degrees of freedom are coupled each other in J^^j. Among 
the several interactions, the inter-site Coulomb repulsion takes 
the largest energy scale. Thus, with lowering temperature {T), 
the charge sector will be frozen, at first. The expected spin or- 
dering temperature is larger than that for the orbital because 
of the large magnitude of the spin operator, 2 and 5/2. 

First we focus on origin of the CO transition accompanied 
by the electric polarization. We pay our attention to the polar 
CO state characterized by the momentum q = (l/3, 1/3,0) = 
q3 [see in Fig. 1(b)]. As suggested by Yamada and cowork- 
ers Jgt], this CO, which we call CO1/3, shows the electric po- 
larization due to the charge imbalance between the triangu- 
lar lattices. Along [110], charge alignments are schematically 
• o • • o» in the upper layer, and • o o • 00 in the lower one, 
where o (•) represents Fe^+ (Fe^+). This CO1/3 competes 
with the non-polar CO with q = (1/2,0,0) =q2 (CO1/2), and 
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FIG. 2: (a): A mean-tield phase diagram at finite T calculated in 
J^. (b): Charge correlation function /V(q) at VcNNN/^abNN = 0-6. 

(c) : A'(q) at VcNNN/VabNN =0.5 (bold lines) and 0.7 (broken lines). 
The insets show schematic views of the CO1/4 and CO1/2 structures. 

(d) : Electric polarization P for several VcNNN/KbNN values. A value 
of VcNN/VabNN is chosen to be 1.2 in (a)-(d). 



that with (1/4, 1/4, 1) = q4 (CO1/4) [see insets of Fig. 2(c)]. 
In both of COi/2 and COi ^4, equal numbers of Fe^+ and Fe^+ 
exist in the lower and upper layers, and the electric polariza- 
tion does not turn up. 

To examine stability of the polar CO 1/3, we analyze 
in a pair of the triangular lattices, which is the minimal unit 
for the electric polarization. We introduce the largest three in- 
teractions Vij [see Fig. 1(a)]: the inter-layer NN interaction 
(Vcnn). the intra-layer NN one (VabNN) and the inter-layer 
next-nearest neighbor (NNN) one (VcNNn)- When the 1/r- 
type Coulomb interaction is assumed, we have VcNN /VabNN = 
1.2 and VcNNN /VabNN = 0.77 for LuFe204. In the case of 
VcNN — VcNNN = 0, the above three CO's are the degenerate 
ground states (GS). Before going to precise calculations, we 
have performed the mean-filed (MF) calculation to obtain a 
global phase diagram. At T —Q, CO 1/2 (CO 1/4) appears in 
the region of 2Vcnnn > Knn (2Vcnnn < Knn), and only on 
the phase boundary, CO1/3 appears. At finite T [Fig. 2(a)], we 
have found that CO 1/3 is stabilized in a wide region between 
COi/2 and CO1/4. 

Unbiased calculations have been performed by utilizing the 
Monte Carlo (MC) simulation. A pair of the triangular lattices 
of L X L sites (L =6, 12) with periodic-boundary condition 
is analyzed by the multi-canonical MC methods |17], where 
8 X 10*' MC steps are used for measurement. The resulting 
charge-correlation functions A^(q) = i^^Li; (6i6/)^"' "^'^ 
[Figs. 2(b) and (c)] are qualitatively consistent with the MF 
phase diagram. At Knnn /VabNN = 0.6 [Fig. 2(b)], A^(q3) 
shows a hump around T/VabNN = 0.18 and keeps the largest 
value down to low T. On the other hand, at Vcnnn /VabNN = 
0.5(0.7) [Fig. 2(c)], a dominant charge correlation function 
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FIG. 3: Charge correlation function A'(q) (a), electric polarization 
P (b), and spin correlation function 5(q) (c) calculated in .Ji^y + -^- 
A schematic view of the charge and spin structure below the ferri- 
magnetic ordering temperature is shown in (d). Spin directions at the 
sites marked by broken circles are not uniquely determined. 



changes from A^(q3) to A^(q2) [A^(q4)] with decreasing T . 
This sequential phase transition may explain the experiments 
in YFe204_5; the observed superlattice spots in the electron 
diffraction change from (n/3 «/3 0) to («/4 n/4 0) with 
decreasing T llisll . The spontaneous electric polarization is 
simulated by calculating P = (p^)'''^ [Fig. 2(d)]. We define 
P = -L"' (L,ec/ - Let) 2' where Y.ieu(L) represents a summa- 
tion of sites in the upper (lower) layer. At VcNNN /VabNN=0.6, a 
humpinP around r/VabNN=0.18 corresponds to that inA^(q3) 
II19II . Apart from VcNNN /VabNN=0.6, a reduction of P is corre- 
lated with the increases of A^(q2) and A'(q4). 

The stability of the polar CO1/3 is caused by the charge 
frustration and an entropy gain on the paired-triangular lat- 
tices. Focus on the Fe sites located on the dotted Unes in 
Fig. 1(b). The charge alignment at the Fe sites is respon- 
sible for the electric polarization. Since these sites are sur- 
rounded by three Fe^+ and three Fe^+ in the plane, the in- 
plane Coulomb interactions VabNN are canceled out. Thus, the 
charge imbalance between the two layers occurs easily with- 
out loss of VabNN, and the polar CO1/3 is stabilized due to 
the inter-layer Coulomb interactions. At these Fe sites, large 
charge fluctuation remains at finite T and contributes to an en- 
tropy gain. This situation in CO 1/3 is in contrast to CO1/2 and 



COi/4 where all sites are equivalent. 

Now we examine coupling between the electric polariza- 
tion and the magnetic ordering. The Hamiltonian ,yfv + 
is analyzed in a pair of the Lx L triangular lattices by utiliz- 
ing the multi-canonical MC method. The spin operators are 
set to be Ising spins, and the exchange constants are estimated 
by using the energy parameters for LaFe03 [20]. The orbital 
PS's in are set to be zero. This is reasonable, because the 



expected SO temperature is much higher than that for orbital 
due to the large magnitudes of the spin operators. The result- 
ing P, N{q) and the spin correlation functions S{q) are plotted 
in Fig. 3. Around r/yabNN=0.18 in Fig. 3(a), A^(q3) shows a 
shoulder, corresponding to a hump in A^(q3) which is calcu- 
lated without J^j. Further increase in A^(q3 ) occurs around the 
spin ordering temperature r/yabNN=0-l [see Fig. 3(c)]. The 
temperature dependence of P [Fig. 3(b)] almost follows that 
of A^(q3). These results clearly show that the SO strongly sta- 
bilizes COi/3 accompanied by the electric polarization. The 
resulting SO is a ferrimagnetic structure characterized by the 
momentum q3 [Fig. 3(d)], which is consistent with the neu- 
tron diffraction experiments 111 ll 11211 . 

The remarkable enhancement of the electric polarization 
below the SO temperature is caused by the spin-charge cou- 
pling and the spin frustration. Focus on the Fe^+ sites in the 
2Fe^+-Fe^+ (upper) layer in Fig. 3(d). Since these Fe ions 
are suiTounded by the NN three-up and three-down spins of 
Fe^+, spin directions at the sites are not determined uniquely. 
There is a large number of degenerate spin states, which con- 
tribute to an entropy gain. Because this spin structure is real- 
ized on the polar CO1/3, this CO accompanied by the electric 
polarization is reinforced through the spin-charge coupling in 
J^j. Essence of this mechanism is based on the spin frustra- 
tion due to the antiferromagnetic configuration for the Fe^+ 
spins, and does not depend on detailed parameter choice. The 
present observation indicates a possibility of the large ME re- 
sponse ifloilTsi lT?]. By applying magnetic field in the CO/SO 
phase, melting of CO 1/3 and reduction of the electric polar- 
ization are expected. 

Finally, we pay our attention to the orbital state. Until now, 
there are no experimental reports for the long-range orbital or- 
der and/or the Jahn-Teller lattice distortion in /?Fe204. Thus, 
here we examine a possible orbital state in the CO/SO struc- 
ture. By using the relation ^/ nf — 0, Hamiltonian Jfj is 
mapped onto an effective orbital model defined on the Fe^+ 
honeycomb lattice in the 2Fe^+-Fe'^+ layer [see Fig. 3(d)]. 
The model is easily derived as 



Jif-c — —J 



i+5„ 
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Here, 5/ indicates a NN Fe-Fe bond labeled by direc- 
tions I = {a,l5,Y) [see Fig. 4(a)], and J{> 0) is the ex- 
change constant. We redefine the orbital PS operator 
as t! cos{K/2 + 2Kni/3)T.'' + sm{K/2 + 2Kni/3)Tf with 
{na,np,ny) — (0, 1,2). This operator represents a projection 
of Tj on the NN Fe-Fe bond direction. When (a,l5,Y) is re- 
placed by {x,y,z) in the Cartesian coordinate, this model cor- 
responds to the e^,-orbital model on the cubic lattice ||21, 22] 
which has been studied for KCUF3 and LaMn03. A related 
model is examined from the view point of the quantum com- 
putation i23ll . It is worth noting that the interaction between 
orbitals explicitly depends on the bond direction. Although 
the honeycomb lattice is bipartite, there are intrinsic frustra- 
tion effects; when PS's are arranged to gain bond energy for 
one direction, these are not fully favorable for other bonds. 
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FIG. 4: (a): A schematic picture of the Fe^+ sublattice in the 2Fe^+- 
Fe-'^ layer. Bold arrows and bold bars represent the orbital PS direc- 
tion and the PS projection components T? on the NN bonds, respec- 
tively, (b): One of the ground state PS configurations, (c): Specific 
heat C{T) for several cluster sizes. Inset shows a magnification for a 
low T region. Broken line represents (cos^ 30) '/^ in L = 9. (d): One 
of the PS configurations at finite temperature below Tq- 



This model provides a non-trivial orbital state. The Hamil- 
tonian is rewritten as ^ = (■//2)L(,7)(t/ - Tj)^ - JLnirff, 
where the second term becomes a numerical constant. This 
form indicates a large number of the classical degenerate GS 
which satisfy t' = rj for each NN ij pair. One of the PS con- 
figurations in the GS is shown in Fig. 4(b). The large de- 
generacy originates from the unusual orbital interaction. The 
momentum representation for the effective orbital interaction 
is given as = £k V^(k)'/(k) v/(k) with a vector \(f{k) = 
[T^{k),T~{k),T^{k),mk)] and a 4 X 4 matrix /(k). Here, A 
and B represent the two sublattices on the honeycomb lattice. 
Surprisingly, the lowest eigen-value for 7(k) is a momentum- 
independent flat band at —37/4. Thus, all eigen states with 
different k provide the degenerate GS. At finite T, we ana- 
lyze J^^t: on a 2 X L X L site cluster (L = 2 ^ 12) by using 
the multi-canonical MC method, where the PS operators are 
treated as classical vectors in the T^-T~ plane. In the specific 
heat C{T) [Fig.4(c)], a peak appears around T /J — 0.025 for 
L = 3. This characteristic temperature, denoted as Tq, is much 
lower than the mean-field ordering temperature 37/8. The 
peak shifts to lower T and becomes sharp with increasing L. 
We also calculate the orbital correlation functions 3^'"' (q) = 
L-^'£ije''i<''--''j'>{T/Tj") {l,m = for all possible q's in 
clusters. However, the magnitudes are less than 20% of their 
maximum values, and remarkable increases in (q) are not 
observed with increasing L. To confirm a precise picture for 
the orbital state below Tq in the thermodynamic limit, numer- 
ical calculations on larger sizes are required. However, within 
the present simulations, we have found a hidden parameter. 



(cos230)i/2 = ([(2L2)-i^.cos30,]2)1/2^ with the PS angle 
di[= t&rT^ [—T^ /Tl')]. This parameter grows up below 7b and 
approaches to 1 [see the inset of Fig. 4(c)]. This result implies 
that 0, at each site takes one of the three angles {2nn)/3 with 
n = (0 , 1 , 2 ) , or one of the three (2 TTn + TT ) /3 . One configura- 
tion is shown in Fig. 4(d) where (2L^)^' cos 30/ = +1. At 
r = 0, any values for 0,- are allowed, as long as the relation 
t' = rj is satisfied. Values of 0; are fixed to be {2Kn)/3 or 
(27i;n + 7r)/3 by thermal fluctuations at finite T. This is the 
so-called order by fluctuation mechanism. 

Summarizing, we present a theory of /?Fe204 as an exotic 
dielectric material with the CO driven electric polarization 
and the ME effects caused by the spin-charge coupling. The 
charge frustration allows the charge imbalance without inver- 
sion symmetry, i.e. the electric polarization. This CO is re- 
inforced by the magnetic ordering where the spin frustration 
contributes to the entropy gain. The present mechanism of 
the ME effects is entirely different from the spin-lattice cou- 
pling mechanism proposed in the previous multiferroics, e.g. 
/JMnOs. A possible orbital state is examined in the CO/SO. 
The effective orbital interaction does not depend on the mo- 
mentum, indicating divergence of fluctuations. Instead of a 
conventional long-range order, a kind of the angle order of the 
pseudo spins grows up. The present study provides an insight 
for searching of a new class of "electronic" multiferroics in 
correlated electron systems. 
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